function [area,circum]=circle(rad)
area=pi.*rad.^2;
circum=2.*pi.*rad;
end